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Abstract 



Results are reported of a search for a deviation in the jet production cross section from 
the prediction of perturbative quantum chromodynamics at next-to-leading order. 
The search is conducted using a 7 TeV proton-proton data sample corresponding to an 
integrated luminosity of 5.0 fb^^, collected with the Compact Muon Solenoid detector 
at the Large Hadron Collider. A deviation could arise from interactions characterized 
by a mass scale A too high to be probed directly at the LHC. Such phenomena can be 
modeled as contact interactions. No evidence of a deviation is found. Using the CLg 
criterion, lower limits are set on A of 9.9 TeV and 14.3 TeV at 95% confidence level for 
models with destructive and constructive interference, respectively. Limits obtained 
with a Bayesian method are also reported. 
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1 Introduction 

Interactions at an energy scale much lower than the mass of the mediating particle can be mod- 
eled by contact interactions (CI) [T-'4] governed by a single mass scale conventionally denoted 
by A. A search for contact interactions is therefore a search for interactions whose detailed 
characteristics become manifest only at higher energies. Contact interactions can affect the 
jet angular distributions as well as the jet transverse momentum (pj) spectra, particularly for 
low-rapidity jets. Lower limits on A have been set by the CDF |j5|, DO HI, and ATLAS |7] collab- 
orations. The Compact Muon Solenoid (CMS) collaboration has previously measured the dijet 
angular distribution [8J using a data set of -^s = 7TeV proton-proton collisions correspond- 
ing to an integrated luminosity of 2.2 fb^^, and found A > 8.4 TeV and A > ll.ZTeV at 95% 
confidence level (CL), for models with destructively and constructively interfering amplitudes, 
respectively. 

The inclusive jet px spectrum, i.e., the spectrum of jets in p + p — )■ jet + X events, where X 
can be any collection of particles, is generally considered to be less sensitive to the presence 
of contact interactions than the jet angular distribution. This perception is due to the jet px 
spectrum's greater dependence on the jet energy scale (JES) and on the parton distribution 
functions (PDF), which are difficult to determine accurately. However, considerable progress 
has been made by the CMS collaboration in understanding the JES [|9J. The understanding of 
PDFs has also improved greatly at high parton momentum fraction ||T0] - IT2| , in part because of 
the important constraints on the gluon PDF provided by measurements at the Tevatron IIT3l[T4| . 
These developments have made the jet pj spectrum a competitive observable to search for 
phenomena described by contact interactions, reprising the method that was used in searches 
by CDF [15J and DO [16J. 

In this paper, we report the results of a search for a deviation in the jet production cross section 
from the next-to-leading-order (NLO) quantum chromodynamics (QCD) prediction of jets pro- 
duced at low-rapidity with transverse momenta >500GeV. The analysis is based on a 7 TeV 
proton-proton data sample corresponding to an integrated luminosity of 5.0 fb^^, collected 
with the CMS detector at the Large Hadron Collider (LHC). 

2 Theoretical models 

The experimental results are interpreted in terms of a CI model described by the effective La- 
grangian t3nl7| 

^ = C^(nLl''^L)mLlF^L), (1) 

where denotes a left-handed quark field and ^ = +1 or — 1 for destructively or construc- 
tively interfering amplitudes, respectively. The amplitude for jet production can be written 
as 

« = «SM + A flci 

where asm arid uqi are the standard model (SM) and contact interaction amplitudes, respec- 
tively. Since the amplitude is linear in A = 1 / A^, the cross section Ck in the kth jet pj bin is 
given by 

a,, = Ck + hA + ak A^, (2) 
where c^, bj^, and a^^ are jet-px-dependent coefficients. 

We use models characterized by the cross section QCD^lo + CI(A), where QCD^lo = ^^k is the 
inclusive jet cross section computed at next-to-leading order, and CI(A) = A + flj. A^ parame- 
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2 Theoretical models 



terizes the deviation of the inclusive jet cross section from the QCD prediction arising from the 
hypothesized contact interactions. The QCD^lq cross section is calculated with version 2.1.0- 
1062 of the fastNLO program with scenario table fnl2332y0.tab [18J using the NLO CTEQ6.6 
PDFs |19|. We do not unfold the observed inclusive jet pj spectrum. Instead, the NLO QCD jet 
Pj spectrum is convolved with the CMS jet response function, where the jet energy resolution 
(JER) CTpj for low-rapidity jets is given by 

= p,^-_ + -^ + c^, (3) 

with n = 5.09, s = 0.512, m = 0.325, c = 0.033, and compared directly with the observed 
spectrum using a likelihood function. Equation ^ is the standard form for the calorimeter 
resolution function, modified to account for a weak pj dependence of the coefficient of the 
ipj^) stochastic term and to model better the resolution of low pj jets by using a negative 
coefficient for the (Pj^) noise term. For brevity, we shall refer to the smeared spectrum as the 
NLO QCD jet pj spectrum. 

The signal term CI (A) is modeled by subtracting the leading-order (LO) QCD jet cross section 
(QCDlo) from the LO jet cross section computed with a contact term. The leading-order jet 
Pj spectra are computed by generating events with and without a CI term using the program 
PYTHIA 6.422, the Z2 underlying event tune flZlEQl, and the same CTEQ PDFs used to calculate 
QCDnlo- The generated events are processed with the full CMS detector simulation program, 
based on Geant4 Il2l1 . Interactions between all quarks are included (Appendix [A) and we 
consider models both with destructive and constructive interference between the QCD and 
CI amplitudes. We note that NLO corrections to the contact interaction model have recently 
become available [22J, and we plan to use these results in future studies. These corrections are 
expected to change the results by less than 5%. 

The jet pj dependence of CI(A) is modeled by fitting the ratio / = [QCDi,jlo + CI(A)]/QCDnlo 
simultaneously to four pythia CI models with A = 3,5,8, and 12 TeV. The fit is performed in 
this manner in order to construct a smooth interpolation over the four cross section ratios. Sev- 
eral functional forms were investigated that gave satisfactory fits, including the ansatz [|23| : 

^ = '+p^ (lo^)'' i^)^^' iwih/T (i^) • 

In a generator-level study, we verified the adequacy of the extrapolation of Eq. Q up to 25 TeV. 
The results of fitting Eq. (j4]| to models with destructive interference are shown in Figure [l] The 
fit shown in Fig. [ljuses the central values of the JES, JER, and PDF parameters and the renor- 
malization (f/j-) and factorization {}if) scales set to j^t = }i{ = jet pj. Models with constructive 
interference are obtained by reversing the sign of the parameter p^. The fit parameters are 
given in Table [l] Figures |2] and |3] show model spectra in the jet pj range 500 < pi < 2000 GeV 
for values of A that are close to the limits reported in this paper. Figure |2] shows that the jet 
production cross section is enhanced at sufficiently high jet pj. However, for interactions that 
interfere destructively, the cross section can decrease relative to the NLO QCD prediction. For 
example, for A = 10 TeV, the QCDi,jlo+CI cross section is lower than the QCD^lo cross section 
for jet Pt < 1.3 TeV. Figure [3] shows the contact interaction signal, CI(A), as a function of jet pj. 
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Figure 1: The cross section ratios, / = [QCD^lo + CI(A)]/QCDi,jlo/ with A = 3,5,8, and 
12 TeV. The points with error bars are the theoretical values of the cross section ratios. The 
curves are the results of a fit of Eq. Q simultaneously to the four cross section ratios. The 
NLO QCD jet px spectrum is calculated using the nominal values of the JES, JER, PDF, renor- 
malization and factorization scales for models with destructive interference. The values of the 
parameters of the fit are given in Table [l] 



Table 1: The fit parameters associated with Fig. [ij The first row lists the values of the parame- 
ters pi, p2/ P3/ and p4, while the remaining rows list the elements of the associated co variance 
matrix. 
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3 Experimental setup 

The CMS coordinate system is right-handed with the origin at the center of the detector, the 
X axis directed toward the center of the LHC ring, the y axis directed upward, and the z axis 
directed along the counterclockwise proton beam. We define (p to be the azimuthal angle, 9 to 
be the polar angle, and the pseudorapidity to be ?/ = — ln[tan(0/2)]. The central feature of 
the CMS apparatus is a superconducting solenoid of 6 m internal diameter, operating with a 
magnetic field strength of 3.8 T. Within the field volume are the silicon pixel and strip trackers 
and the barrel and endcap calorimeters with \ri\ < 3. Outside the field volume, in the forward 
region, there is an iron/ quartz-fiber hadron calorimeter (3 < \ri\ < 5). Further details about 
the CMS detector maybe found elsewhere ||24J. 

Jets are built from the five types of reconstructed particles: photons, neutral hadrons, charged 
hadrons, muons, and electrons, using the CMS particle-flow reconstruction method |25| and 
the anti-A:x algorithm with a distance parameter of 0.7 ||26l - |28| . The jet energy scale correction 
is derived as a function of the jet px and rj, using a px-balancing technique ID, and applied to 
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Figure 2: The cross section ratios, / = [QCDnlo + CI(A)]/QCDnlo/ with A = 8, 10, 12, and 
14 TeV, for models with destructive (left) and constructive (right) interference. 
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Figure 3: The CI signal spectra, defined as daQci:)+ci/dpj — dcTQCD/dfj (pb/GeV) with A = 
8, 10, 12 and 14 TeV, for models with destructive (left) and constructive (right) interference. 



all components of the jet four momentum. 

The results reported are based on data collected using un-prescaled single-jet triggers with pj 
thresholds that were changed in steps from 240 to 300 GeV during the data-taking period. The 
trigger thresholds were changed in response to the increase in instantaneous luminosity. The 
jet trigger efficiency is constant, ~98.8%, above ~400 GeV, well below the search region. Events 
with hadron calorimeter noise are removed |29| and each selected event must have a primary 
vertex within 24 cm of the geometric center of the detector along the z axis and within 0.2 cm 
in the transverse x-y plane, defined by criteria described in 1130 1 . The search is restricted to 
\f]\ < 0.5 where the effects of contact interactions are predicted to be the largest llH4l. The jet 
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Pt spectrum is divided into 20 pj bins in the search region 507 < pj < 2116 GeV, where the 
bin width is approximately equal to the jet resolution (Tpj given in Eq. (|3|. No jets are observed 
above 2000 GeV transverse energy. 



4 Results 



In Figure |4] we compare the observed inclusive jet pj spectrum with the NLO QCD jet pj spec- 
trum, which is normalized to the total observed jet count in the search region using the normal- 
ization factor 4.007 ± 0.009 (stat.) fb"^ (Section|5). The normalization is the ratio of the observed 
jet count to the predicted cross section in the search region. The data and the prediction are in 
good agreement as indicated by two standard criteria, the Kolmogorov-Smimov probability 
Pr(KS) of 0.66 and the per number of degrees of freedom (NDF) of 23.5/19. Table [2] lists 
the observed jet counts. Figure |5] compares the observed jet pj spectrum in the search region 
with model spectra for different values of A, for models with destructive interference. Figure |6] 
compares the data with models with constructive interference. 
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Figure 4: The observed jet pj spectrum compared with the NLO QCD jet pj spectrum (left). The 
bands represent the total uncertainty in the prediction and incorporate the uncertainties in the 
PDFs, jet energy scale, jet energy resolution, the renormalization and factorization scales, and 
the modeling of the jet pj dependence of the parameters in Eq. (Q. The ratio of the observed 
to the predicted spectrum (right). The error bars represent the statistical uncertainties in the 
expected bin count. 



5 Statistical analysis 

Since there are no significant deviations between the observed and predicted spectra, the re- 
sults are interpreted in terms of lower limits on the CI scale A using the models described in 
Section |2] The dominant sources of systematic uncertainties are associated with the JES, the 
PDFs, the JER, the renormalization (^r) arid factorization (f/f) scales, and the modeling param- 
eters of Eq. Non-perturbative corrections are less than 1% for transverse momenta above 
~400 GeV [30J, are negligible compared with other uncertainties, and are therefore not applied 
to our analysis. 
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5 Statistical analysis 



Table 2: The observed jet count for each jet pj bin in the range 507-2116 GeV. 
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Figure 5: The data compared with model spectra for different values of A for models with 
destructive interference (left). The ratio of these spectra to the NLO QCD jet pi spectrum 
(right). 




In the search region, the inclusive jet spectrum has a range of five orders of magnitude, which 
causes the limits on A to be sensitive to the choice of the normalization factor and the size of the 
data sets. We have found that a few percent change in the normalization factor can cause limits 
to change by as much as 50%. Therefore, for the purpose of computing limits, we have chosen 
to sidestep the issue of normalization by considering only the shape of the jet pi spectrum. This 
we achieve by using a multinomial distribution, which is the probability to observe K counts, 
Ny, / = 1, • • • , K, given the observation of a total count N = Y^=\ ^j- The likelihood is then 
defined by 

''<°|'^-"' = ^^n(^) - (5) 

where K = 20 is the number of bins in the search region, Nj is the jet count in the jth jet px 
bin, D = Ni, • • • , Nk, (t = Y^f^i dj and N are the total cross section and total observed count, 
respectively, in the search region, and the symbol co denotes the nuisance parameters pi, • • ■ , p4 
inEq. Q. 
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Figure 6: The data compared with model spectra for different values of A for models with 
constructive interference (left). The ratio of these spectra to the NLO QCD jet pj spectrum 
(right). 



We account for systematic uncertainties by integrating the likelihood with respect to a nuisance 
prior tt{co). In practice, the likelihood is averaged over the nuisance parameters, co, using a dis- 
crete representation of the prior 7i{co) constructed as described in Section 5.1 This calculation 
yields the marginal likelihood p(D|A) Ylm=i p(D|A,a;m), where M is the number of points 
sampled from the nuisance prior tt{cv) described in Appendix B.l which is the basis of the 
limit calculations. The likelihood functions for models with destructive and constructive inter- 
ference are shown in Figure [7[ 
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Figure 7: The likelihood functions assuming a model with either destructive (left) or construc- 
tive (right) interference. The dashed curve is the likelihood function including statistical uncer- 
tainties only and the central values of all nuisance parameters. The solid curve is the likelihood 
marginalized over all systematic uncertainties. 
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6 Summary 



5.1 Uncertainties 

In principle, a discrete representation of the nuisance prior n{co) can be constructed by sam- 
pling simultaneously the JES, JER, PDFs, and the three values of f/f and fir- Pt/2, pj, and 
Ipj. However, the CTEQ collaboration [19] does not provide a sampling of PDFs. Instead, 
CTEQ6.6 contains 44 PDF sets in which the 22 PDF parameters are shifted by approximately 
±1.64 standard deviations. If we assume the Gaussian approximation to be valid, we can con- 
struct approximate 20 x 20 covariance matrices for the jet spectra from the 44 PDF sets. Us- 
ing these matrices, we generate ensembles of six correlated spectra: QCD^lq/ QCDlq/ ^^id 
(QCD + CI(A))lo with A = 3,5,8, and 12TeV. The generation is performed for models both 
with destructive and constructive interference. The details of our procedure, which also in- 



cludes simultaneous sampling of the JES and JER parameters, are given in Appendix B.l 

For a given set of values for the JES, JER, PDF, fi^, and fif parameters, Eq. Q is fitted to the ratio 
(QCD^LO + CI) / QCD^LO simultaneously to the four models with A = 3, 5, 8, and 12 TeV. We 
then sample a single set of the four nuisance parameters co = pi, p2, ps, P4 from a multivariate 
Gaussian using the fitted values and the associated 4x4 covariance matrix. The sampling 
and fitting procedure is repeated 500 times, thereby generating a discrete representation of the 
nuisance prior rc{co) that incorporates all uncertainties. We have verified that our conclusions 
are robust with respect to variations in the size of the sample that represents tt{<x>). 

5.2 Lower limits on A 

We use the CLg criterion ||3T1 132| to compute upper limits on A. For completeness, we give 



the details of these calculations in Appendix B.2 Using the procedure described in the Ap 



pendix, we obtain 95% lower limits on A of 9.9 TeV and 14.3 TeV for models with destructive 
and constructive interference, respectively. These more stringent limits supersede those pub- 
lished by CMS based on a measurement of the dijet angular distribution ||8l. The current search 
is more sensitive than the earlier dijet search as evidenced by the expected limits, which for this 
analysis are 9.5 ± 0.6 TeV and 13.6 ±1.6 TeV, respectively, obtained using 5 fb^^ of data. 

Limits are also computed with a Bayesian method (Appendix |B.3| using the marginal likeli- 
hood p(D|A) and two different priors for A: a prior flat in A and a reference prior [,,33'- 35] . 
Using a flat prior, we find lower limits on A of 10.6 TeV and 14.6 TeV at 95% confidence level for 
models with destructive and constructive interference, respectively. The corresponding limits 
using the reference prior are 10.1 TeV and 14.1 TeV at 95% confidence level, respectively. 



6 Summary 

The inclusive jet pj spectrum of 7 TeV proton-proton collision events in the ranges 507 < px < 
2116 GeV and |//| < 0.5 has been studied using a data set corresponding to an integrated lu- 
minosity of 5.0 fb^^. The observed jet px spectrum is found to be in agreement with the jet 
Px spectrum predicted using perturbative QCD at NLO when the predicted spectrum is con- 
volved with the CMS jet response function and normalized to the observed spectrum in the 
search region. Should additional interactions exist that can be modeled as contact interactions 
with either destructive or constructive interference, their scale A is above 9.9 TeV and 14.3 TeV, 
respectively, at 95% confidence level. We plan to extend this study to the full 8 TeV CMS data 
set, making use of a recently released program [36J to calculate at next-to-leading order the 
inclusive jet pj spectrum with contact interactions. 

It is noteworthy that the limits reported in this paper, which are the most sensitive limits pub- 
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lished to date, have been obtained reprising the classic method to search for contact interac- 
tions: namely, searching for deviations from QCD at high jet transverse momentum. 
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A PYTHIA 6.422 contact interaction configuration 

The scale A is defined by the CI model in pythia. In order to facilitate the re-interpretation 
of the results using a different model, we provide the details of the pythia configuration in 
Table |3] for A = 8 TeV and final state parton transverse momenta, pj, in the range 170 < pj < 
230 GeV. 
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B Statistical details 



Table 3: pythia 6.422 configuration for A = 8 TeV contact interactions. 



PYTHIA 6.422 settings specific to contact interactions 



Settings 






Description 


ITCJVl(5)=2 


Switch on 


contact int. for all quarks 


RTCM(41)=8000 




Set contact scale A to 8 TeV 


RTCM(42)=1 




Sign of contact int. is + 


MSUB(381)=1 




q/qy 


via QCD plus a contact int. 


MSUB(382)=1 




mk 


via QCD plus a contact int. 


MSUB(13)=1 




q^q/ - 


gg via normal QCD 


MSUB(28)=1 




q/g^ 


► qjg via normal QCD 


MSUB(53)=1 




gg^ 


q^tq^^ via normal QCD 


MSUB(68)=1 




gg^ 


► gg via normal QCD 


CKIN(3)=170 




minimum pj for hard int. 


CKIN(4)=230 




maximum pj for hard int. 



B Statistical details 

B.1 Nuisance prior 

We approximate the nuisance prior t[{co) starting with two sets of ensembles. In the first, the 
six 20-bin model spectra QCDnlo/ QCDlq, and [QCD + CI(A)]lo with A = 3, 5, 8, and 12 TeV 
are varied, reflecting random variations in the PDF parameters as well as random choices of the 
three fir and ^{ scales, while keeping the JES and JER parameters fixed at their central values; 
we call these the PDF ensembles. In the second set of ensembles, the JES and JER parameters 
are varied simultaneously, while keeping the PDF parameters fixed to their central values and 
the renormalization and factorization scales at their nominal values; we call these the JES/ JER 
ensembles. 

B.I .1 Generating the PDF ensembles 

In the PDF ensembles, each of the six model spectra is sampled from a multivariate Gaussian 
distribution using the associated 20 x 20 covariance matrix. For each model spectrum, the 
covariance matrix is approximated by 

22 22 

Cnm = ^-^"i ^^mj, (6) 

/=1;=1 

where AX„, = (X+ — X^^)/2 and X^- are the cross section values for nth jet bin associated with 
the + and — variations of the fth pair of CTEQ6.6 PDF sets. CTEQ [19J publishes approximate 
90% intervals. We therefore approximate 68% intervals by dividing each AX by 1.64. The 
correlation induced by the PDF uncertainties across all six model spectra is maintained by 
using the same set of underlying Gaussian variates during the sampling of the spectra. 

B.I .2 Generating the JES/JER ensembles 

In the JES/JER ensembles, the JES and JER parameters are sampled simultaneously for the five 
model spectra QCDlq/ and (QCD + CI)lo with A = 3,5,8, and 12 TeV, yielding ensembles 
of correlated shifts from the central JES, JER, and PDF values of the QCDlq and (QCD+CI)lo 
spectra. For example, we compute the spectral residuals 5a = QCD' — QCD^-gj^t^ai/ where QCD' 
is the shifted jet pj spectrum and QCD^g^tj-ai is the jet px spectrum computed using the central 
values of the JES, JER, and PDF parameters. Coherent shifts of the jet energy scale are calculated 



B.2 CLs calculation 
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for every jet in every simulated event. The jet pj is shifted by x3 for each component of the 
jet energy scale uncertainty, of which there are sixteen, where x is a Gaussian variate of zero 
mean and unit variance, and ^ is a jet-dependent uncertainty for a given component. The 
contributions from all uncertainty components are summed to obtain an overall shift in the jet 
Pj. From studies of dijet asymmetry and photon+jet pj balancing, the uncertainty in the jet 
energy resolution is estimated to be 10% in the pseudorapidity range < 0.5 f30]. We sample 
the jet energy resolution using a procedure identical to that used to sample the jet energy scale, 
but using a single Gaussian variate. 

B.I .3 Generating the JES/JER/PDF ensemble 

Another ensemble is created, from the PDF ensembles and the JES/JER ensembles, that ap- 
proximates simultaneous sampling from the JES, JER, PDF, renormalization, and factorization 
parameters. We pick at random a correlated set of six spectra from the PDF ensembles, and 
a correlated set of five spectral residuals from the JES/JER ensembles. The JES/JER spec- 
tral residuals Scr are added to the corresponding shifted spectrum from the PDF ensembles, 
thereby creating a spectrum in which the JES, JER, PDF, ji^, and }i{ parameters have been ran- 
domly shifted. The NLO QCD spectrum (from the PDF ensembles) is shifted using the LO QCD 
JES/JER spectral residuals in order to approximate the effect of the JES and JER uncertainties 
in this spectrum. 

The result of the above procedure is an ensemble of sets of properly correlated spectra QCD^lo + 
CI(A) with A = 3,5,8, and 12TeV, in which the JES, JER, PDF, ji^ and f/f parameters vary 
randomly. The ansatz in Eq. ^ is then fitted to the quartet of ratios [QCD^lo + CI(A)] / 



QCDi,jLo as described in Section 5.1 to obtain parameter values for pi, p2, ps, and p4. Five hun- 
dred sets of these parameters are generated, constituting a discrete approximation to the prior 

n{w) = 7r(pi,p2,p3,P4)- 

B.2 CLs calculation 

Since CLg is a criterion rather than a method, it is necessary to document exactly how a CLg 
limit is calculated. Such a calculation requires two elements: a test statistic Q that depends on 
the quantity of interest and its sampling distribution for two different hypotheses, here A > 0, 
which we denote by H^, and A = 0, which we denote by Hq. H\ is the signal plus background 
hypothesis while Hq is the background-only hypothesis. For this study, we use the statistic 

Q(A) = t{D,\) = -21n[p(D|A)/p(D|0)], (7) 

where p(D|A) is the marginal likelihood 

p(D|A) = j p(D|A,a;) 7r(a;) da;, 

1 M (8) 
~ a7 E p(D|A,a;„0, 

m=l 

where M = 500 is the number of points on = pi, p2, ps, P4 sampled from the nuisance prior 
n{co) described in Appendix |B.1[ We compute the sampling distributions 



and 



P(Q|Ha) = j 5[Q-t{D,\)]p{D\\)dD, (9) 
p(Q|Ho) = |^[Q-f(D,A)]p(D|0)dD, (10) 
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pertaining to the hypotheses Hx and Hq, respectively, and solve 



CLs = p(A)/p(0) =0.05, 



(11) 



to obtain a 95 % confidence level (CLg) upper limit on A, where the p-value p(A) is defined by 



and Qo is the observed value of Q. 

In practice, the CLg limits are approximated as follows: 

1. Choose a value of A, say A*, and compute the observed value of Q, Qo(A*). 

2. Choose at random one of the M = 500 sets of nuisance parameters pi, p2, ps, and p4. 

3. Generate a spectrum of K = 20 counts, D, according to the multinomial distribution, 
Eq. with A = A*, which corresponds to the hypothesis H^. Compute Q = t{D, A*) and 
keep track of how often Q(A*) > Qo(A*). Call this count n\. 

4. Generate another set of 20 counts, D, but with A = 0, corresponding to the hypothesis Hq. 
Compute Q = t{D,\*) and keep track of how often Q(A*) > Qo(A*). Call this count tiq. 

5. Repeat 25, 000 times steps 2 to 4, compute CLg ~ / no and report A = A* as the upper 
limit on A at 95% CL if CLg is sufficiently close to 0.05; otherwise, keep repeating steps 1 
to 4 with different values of A. The algorithm starts with two values of A that are likely 
to bracket the solution and the solution is found using a binary search, which typically 
requires about 10 to 15 iterations. 

B.3 Bayesian calculation 

The Bayesian limit calculations use the marginal likelihood, Eq. (j8|, and two different (formal) 
priors 7r(A): a prior flat in A and a reference prior ||33ti35l , which we calculate numerically II35I . 
An upper limit on A, A*, is computed by solving 



p(A) = Pr[Q(A) > Qo(A)], 



(12) 




(13) 



where p{D) is a normalization constant. The integrals are performed using numerical quadra- 
ture. 
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